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ABSTRACT 

The isolation, simple apparent structure, and low luminosity of the nearby spiral 
galaxy NGC 6503 make it an ideal candidate for an in-depth kinematic and photo- 
metric study. We introduce a new publicly available code, DiskFit, that implements 
procedures for fitting non-axisymmctrics in cither kinematic or photometric data. We 
use DiskFit to analyze new Ha and CO velocity field data as well as Hi kinemat- 
ics from Greisen et al. to search for non-circular motions in the disc of NGC 6503. 
We find NGC 6503 to have remarkably regular gas kinematics that are well-described 
by rotation only. We also use DiskFit and a new l^s-band image of NGC 6503 to 
constrain photometric models of the disc, bar and bulge. We find the galaxy to be 
photometrically dominated by the disc. We find NGC 6503 to contain a faint bar and 
an exponential bulge which together contribute only ~ 5% of the total galaxy light. 
The combination of our kinematic and photometric DiskFit models suggest that 
NGC 6503 contains a weak, end-on bar that may have produced its Type II surface 
brightness profile but is unlikely to be responsible for its strong er-drop. 

Key words: galaxies: kinematics and dynamics — galaxies: structure 



1 INTRODUCTION 

NGC 6503 is a low- luminosity, late-type spiral galaxy that 
is relatively isolated from systems of similar mass. It has 
been the focus of several recent studies to characterize its 
underlying d isc structure. Using Very Large Array (V L.A[]]) 
observations, iGreisen. Spekkens. fe van Moorsel |2009l ) find 
the kinematics of the Hi disc to be remarkably regular. 
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They also find that the d ata require a disc with both 
thin a nd thick components. iPuglielli. Widrow. fc Courteaul 
|2010T l construct dynamical models of NGC 6503 using 
Bayesian statistics and Markov chain Monte Carlo (MCMC) 
techniques. Their extensive modeling places interesting con- 
straints on the structure of the galaxy disc and bulge, mass- 
to- light ratios, and the dark matter halo. They find the bulge 
of the galaxy to have a low Sersic index and a mass-to-light 
ratio lower than the disc, and that the dark matter halo is 
cuspy. 



iFreeland et~aT] (|2010l . hereafter F10) present WIYN0 
High-resolution Infrared Camera (WHIRC) i/-band 
(1.6/un) imaging, in combination with multi- wavelength 
archival imaging and data from the literature, to investigate 
the origin(s) of the observed star-forming ring and central 
a-drop. They find a plateau in the surface brightness profile 



2 Based on observations obtained at the WIYN Observatory. 
The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatory. 
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Figure 1. Observed velocity fields of NGC 6503. All colourscales ar e identical, and all velocity fields have been rotated 30° to the west. 
Note the spatial extent is different in each panel. 1" 0.023 kpc jMould et al.ll2000h . Top left: DensePak Ho data. The diameter of 
the DensePak fibres is 3". Top right: SparsePak Ha data. The SparsePak fibres are 5" in diameter. Bottom left: CO data. The spatial 
resolution is 4.2". Note the smaller spatial e xtent of the CO data compared to the Ho data. Bottom right: Hi vel ocity field derived using 
the Modified Envelope Tracing technique of iGentile et al.l l l2004h on the Hi data cube from [Greisen et al.l | |2009| ). The spatial resolution 
is 14". Note that the spatial extent of the Hi velocity field is about 7 times larger than the optical Ha data. (A colour version of this 
figure is available in the online journal.) 



(Freeman Type II behaviour, iFreemar] Il970h and several 
structures in the disc in addition to the star-forming ring: 
an end-on bar, a circumnuclear di sc, and a nuclear spiral. 
Comp aring to the simulations of iBureau &i Athanassoulal 
(|2005l . hereafter B05), they argue that the bar in NGC 6503 
is strong and that it is this structure which can explain 
the photometric and kinematic features. In this model, 
the remarkable regularity of the kinematics in NGC 6503 
therefore belies a strong non-axisymmetry near the centre. 

Because of its suitability for detailed structural studies, 
we have acquired new, high-quality, multi-wavelength kine- 
matics and photometry of NGC 6503. In this paper, we use a 
newly released code, DiskFit, to perform an in-depth study 
of the underlying structure of the galaxy. We investigate the 
gas kinematics and search for non-circular motions in Ha, 
CO, and Hi velocity field data. We also use A's-band imag- 
ing to constrain photometric models of the disc, bar, and 
bulge. 

In Section [2] we present the new observational data for 
NGC 6503. In Section [3] we introduce and describe the first 
DiskFit release, that we use to fit non-parametric models 
to the galaxy data. The results of the kinematic modeling 
are presented in Section 3] photometric modeling results are 
in Section [5] In Section [(J we discuss the kinematic and 
photometric constraints on the structure of NGC 6503. We 
summarize our results in Section [7] 



2 DATA 

The isolation, simple apparent structure, and low-luminosity 
of NGC 6503 make it an ideal candidate for kinematic and 
photometric study. We therefore compile a multi-wavelength 
dataset of high-resolution spectroscopy and photometry for 
the galaxy. 

We obtain new integral field unit (IFU) Ha veloc- 
ity fields using both the DensePak and SparsePak IFUs 



on the 3.5-m WIYN telescope at Kitt Peak National Ob- 
servatory (KPNO). We also obtain CO(J=1-0) observa- 
tions using the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA0). To minimize the effects of 
beam smearing, we derive a new Hi velocity field from the 
I Greisen et al.l (|2009l ) data cube. We also obtain new Wide- 
field Infrared Camera (WIRCAM) Tfs-band imaging using 
the 3.6-m Canada- France-Hawaii Telescope (CFHlQ). In the 
following sections, we present each of the datasets and de- 
scribe the observations in detail. 



2.1 Ha velocity fields 

We have obtaine d Ha velocity fields of NGC 6503 us- 
ing the Den sePak dBarden. Sawyer, fc Honevcuttll 19981 ) and 
SparsePak jBershadv et al.ll2004 ) IFUs on the WIYN tele- 
scope at KPNO. The DensePak IFU is comprised of 85 work- 
ing 3" diameter fibres arranged in a fixed 43" x 28" rect- 
angle. Similarly, the SparsePak IFU is comprised of 82 5" 
diameter fibres arranged in a fixed 70" x 70" rectangle. 

The DensePak observations were obtained on 2004 July 
30. The instrument setup matched that of ISimon et al.l 



3 Support for CARMA construction was derived from the Gor- 
don and Betty Moore Foundation, the Kenneth T. and Eileen L. 
Norris Foundation, the James S. McDonnell Foundation, the As- 
sociates of the California Institute of Technology, the University 
of Chicago, the states of California, Illinois, and Maryland, and 
the National Science Foundation. Ongoing CARMA development 
and operations are supported by the National Science Founda- 
tion under a cooperative agreement, and by the CARMA partner 
universities. 

4 Based on observations obtained with WIRCam, a joint project 
of CFHT, Taiwan, Korea, Canada, France, and the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the Na- 
tional Research Council (NRC) of Canada, the Institute National 
des Sciences de l'Univers of the Centre National de la Recherche 
Scientifiquc of France, and the University of Hawaii. 
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l|2003f l and ISimon et all (120051 ) . with the Bench Spectro- 
graph in echelle mode with 13 km s" 1 velocity resolution 
over a narrow spectral range centred on Ha. The galaxy was 
observed with 11 DensePak pointings, seven along the major 
axis, two on the minor axis, and one each to the northwest 
and southeast, providing 654 independent velocity measure- 
ments across the disc. Because the Ha emission was quite 
bright, most of the pointings received a single 1200 s expo- 
sure, with a second matching exposure at a few positions 

where the emission was fainter. 

Data r eduction fol l owed the outline of I Simon et al.l 

|2003h and Tsimon et al.l (|2005T l. The data were reduced in 
IRAFjj using the hydra package. The individual frames 
were bias-subtracted and flat fielded, and then cosmic rays 
were removed. The spectrum of each fibre was extracted 
and wavelength calibrated by comparison with a ThAr lamp 
frame. Next, we averaged together the four sky fibres, leav- 
ing out any sky spectra that were contaminated by emission 
lines from NGC 6503. We then removed a linear baseline, 
performed a Gaussian fit to the averaged sky emission near 
Ha, and subtracted the fit from all of the data fibres. Fi- 
nally, individual exposures of the same field were averaged 
together when present, and velocities were measured for each 
fibre by fitting a Gaussian to the observed Ha emission line. 
The median velocity uncertainty from the Gaussian fits was 
1.2 km s . The observed Ha intensities in each pointing 
were cross-corr elated against a narrow -band Ha image of 
the galaxy from lKennicutt et al.l (|2008T ) to determine the ab - 
solute position of each observation fsee lSimon et al.l 12003? ). 
The observed DensePak Ha velocity field is shown in Fig. [1] 
A map of the Ha intensity in the DensePak f ibres is shown 
in Fig. [2] The star-forming ring discussed by iFld is visible 
in the Ha intensity map. 

NGC 6503 was observed with the SparsePak IFU on 
2009 May 16. The STA1 CCD was used with the 316@63.4 
grating in eighth order, centred near Ha, giving a 40 kms -1 
velocity resolution. The SparsePak array was aligned with 
the major axis of the galaxy and three pointings were used 
to cover the length of the galaxy. Individual exposures were 
1200 s, and two exposures were taken at each pointing. A 
ThAr lamp was observed to provide wavelength calibration. 

The SparsePak data were reduced in IRAF using 
the hydra package follow ing the procedure outlined in 
iKuzio de Narav et al.l (|2006l ). The data were bias-subtracted 
and flattened, and the IRAF task dohydra was used to ex- 
tract the spectra. The spectra were wavelength-calibrated 
using a wavelength solution created from the observations 
of the ThAr lamp. The two exposures per pointing were 
combined to remove cosmic rays and to increase the signal- 
to-noise. Sky subtraction was not performed, as the night- 
sky emission lines were used as the reference wavelengths 
(|Osterbrock et akll 19961 ) by which the velocities of the galac- 
tic emission lines were measured. 

The observed fibre velocities were measured by fitting 
Gaussians to both the sky lines and the five galactic emis- 
sion lines of interest: Ha, [Nll]A6548, [Nll]A6584, [Sll]A6717, 



5 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under agreement with the 
National Science Foundation. 
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Figure 2. Ha intensity map from the DensePak observations. 
The dashed whit e ellip se indicates the position of the star- forming 
ring discussed bv lFlOl and in Scction l5.2l The linear flux scaling is 
arbitrary. (A colour version of this figure is available in the online 
journal.) 
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Figure 3. CARMA CO (J=l-0) intensity map. Contour levels are 
at evenly spaced intervals from 10% to 90% of the peak intensity. 
The blue dashed ellip se indicates the position of the star-forming 
ring discussed bv lFldl and in Section l5.2l The beam size is plotted 
in the top left corner. (A colour version of this figure is available 
in the online journal.) 

and [S n] A6731. There was less scatter between the measured 
galactic emission line velocities when using the night-sky cal- 
ibration than when using the ThAr calibration. The arith- 
metic mean of the measured emission-line velocities in each 
fibre was used as the fibre velocity and the error on the fi- 
bre velocity was set to the maximum difference between the 
measured velocities and the mean. Most of the errors are less 
than 5 kms -1 , though a few were greater than ~10 kms . 
The observed SparsePak Ha velocity field is shown in Fig. [T] 

2.2 CO velocity field 

NGC 6503 was observed with the 15-antenna (10-m and 
6-m dishes) CARMA array in March/ April 2008 and Oc- 
tober/November 2009. The instrument setup, data reduc- 
tion and ca l ibrati on are the same as those described in 
iKoda et alj l|201ll ). We used 1642+689 or 1849+670 as 
phase calibrator, interrupting source integrations every 20 
minutes. Absolute flux calibration was done using MARS, 
NEPTUNE or MWC349, whichever was available, and 
3C273/3C345 as bandpass calibrator. We mosaicked the 
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galaxy with 9 pointings: 5 beams along the major axis, and 
2 beams on either side of t he minor axis. Standard routines 
from the miriad package ()Sault. Teuben. fc Wrightl fl995l ) 
were used for calibration and mapping. A final CO datacube 
was produced with a spatial resolution of 4.2" x 3.4", and 
in velocity 5 kms" 1 . 

The CO velocity field shown in Fig.[T]is derived from the 
first moment of the datacube, where only emission exceeding 
2.5 mjy/beam after the data were spatially smoothed by 3" 
and spectrally smoothed by 15 kms -1 were retained. Note 
the smaller spatial extent of the CO data compared to the 
Ha data. In Fig. Owe show the intensity map of the CO data 
an d ind icate the position of the star-forming ring discussed 
bvlFlol. 



2.3 Hi velocity field 

The reduced Hi datacube for NGC 6503 is from 
iGreisen et al.l (|2009l ). Briefly, these data were obtained us- 
ing the VLA in C configuration. There were 127 spectral 
channels, each being separated by 24.4 kHz (5.15 kms -1 ), 
and centred on 26 kms -1 heliocentric radial velocity. The 
total on-source integration time was just over 500 minutes, 
and the data were imaged to produce an angular resolution 
of 14". For complete details on the observing setup and d ata 
reduction, the reader is referred to IGreisen et al.l (|2009T ). 

In order to minimize the effects of beam smearing in 
this high inclination system (see Ta ble [0, we use the M od- 
ified Envelope Tracing technique of Gentile et al, 1 (120041 ') to 
derive a velocity field from the IGreisen et al.l (|2009T ) data 
cube. The resulting velocity field is shown in Fig. [1] Note 
that the radial extent of the Hi data is about 7 times larger 
than that of the Ha data. Beyond r — 300" where beam 
smearing is negligible, the velocity field is almost identical 
to that obtained from the first momen t of the Hi distribu- 
tion in figure 8 of IGreisen et all (|2009r i. At r < 300" where 
beam smearing produces skewed line-of-sight Hi profiles, the 
projected velocities in Fig. [1] have larger amplitudes than in 
the first moment map, by an average of 10 kms - along the 
major axis. 



2.4 Ks-band image 

ifs-band imaging of NGC 6503 was obtained using 
WIRCAM on the CFHT in March 2009. The observations 
were carried out using a 16-point dither pattern repeated on 
each of the four WIRCAM chips, with 3x15 s exposures at 
each position. This resulted in a total of 3 X 16 x 4 = 192 x 
15 s exposures, for a total integration time of 48 minutes. 

Processing of these data, including bad pixel rejection, 
sky subtraction and astrometric and photometric calibra- 
tion, was performed using version 1.9 of the Tiwi pipelin^f). 
The photometric calibration, which is provided by the Tiwi 
pipeline, was manually checked using background 2MASS 
stars in order to confirm its validity. 

We estimated the sky level in the fully-reduced final 
image by calculating the mode of the sky intensities per 
pixel in 20 regions away from the galaxy that were free of 
any other sources. The average and standard deviation of 
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Figure 4. CFHT WIRCAM Ks-band image of NGC 6503, plot- 
ted on a logarithmic scale. The extent of the image is the full 
length of the surface brightness profile discussed in Section I5.3I 
See Fie. UOI for a sky-subtracted, masked view of the inner disc. 
(A colour version of this figure is available in the online journal.) 



the sky values were also computed. The sky fluctuation is 
~ 0.005%, and we therefore subtracted a constant sky level 
across the entire image. The WIRCAM ifs-band image of 
NGC 6503 is shown in Fig. H 



3 DiskFit 

All of the model fits to the data presented in this paper 
were carried out using DiskFit, a newly-released, publicly- 
availablfjj] code that fits non-parametric models to either im- 
ages or velocity fields. While DiskFit will fit axisymmetric 
models, its main purpose is to fit non-axisymmetric mod- 
els either to images or to velocity f ields of disc galaxies, as 
origi nally described by Reese et al.l (|2007 ) for images and by 
both ISpekkens fc Sellwoodl (|200tT ) and lSellwood fc Sanchej 



(|2010l ) for velocity fields. DiskFit supercedes the velfi 
code to model velocity fields, and it represents the first 
public release of the n on-parametri c imag e decomposition 
technique described in iReese et all (|2007h . Both the pho- 
tometric and kinematic branches of DiskFit employ the 
same basic minimization technique, originally described by 
iBarnes fc Sellwoodl (|2003T ). Several new features have been 
added to the algorithm since the original publications: for 
example, DiskFit can now fit for a symmetric outer veloc- 
ity field warp, correct for minor distortions due to seeing, 
handle FITS or text file velocity field inputs, and sparsely 
sample input FITS files. 

As described in detail in previous publica- 



tions 



described m de tail m previous puDlica- 

dReese et al.l I2007I : ISpekkens fc Sellwoodf l2007i : 
ianchezl l 



Scll wood fc Sanche 2010h , DiskFit has several advantages 
over traditional approaches for extracting the kinematic 
or photometric structure of a disc galaxy. The kinematic 
branch of the code, velfit, differs fundam entally from 
the frequently-used ROTCUR (|Begemanlll989T ) "tilted ring" 
algorithm since it fits a specific, physically motivated model 
rather than parametrizing concentric rings of the velocity 
field. By the same token, the photometric branch of the 
code differs fundamentally from popular algorithms such 



7 http : //www. physics . rutgers . edu/~spekkens/diskf it/ 



http : //cf ht .hawaii . edu/Instruments/Imaging/WIRCam/IiwiVersioifllliilxtphVijflf ™. physics . rutgers . edu/~spekkens/velf it/ 
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as Galfit (|Peng et all l2002i . |201Cft in that it fits non- 
parametric disc and bar light profiles rather than specified 
functional forms. DiskFit therefore enables the user to 
model both kinematic and photometric data for nearby 
galaxies with as few assumptions about the parameters of 
the disc and bar as can reasonably be made. Because the 
underlying algorithm is the same for both the kinematic and 
photometric branches of the code, DiskFit is a particularly 
powerful tool for understanding the physical structure of 
galaxies for which both types of data are available. 

DiskFit excels at finding relatively weak, coherent 
asymmetries: it provides a robust estimate of the underly- 
ing disc properties in their presence, and quantitative upper 
limits in their absence. It is also superior to other algorithms 
in that it is capable of providing statistically valid, and real- 
istic, estimates of the uncertainties on the returned param- 
eters using bootstrap realizations of the best-fitting model 
to the data. 

Extensive documentation about the DiskFit approach, 
models a nd capabiliti e s can be found in previou s pub- 
lications jReese et all 120071 ; ISpekkens fc Sellwoodl 120071 ; 
ISellwood fe Sanchej|2O10f i" and on the code website. Below, 
we summarize the salient features of the kinematic and pho- 
tometric models that are relevant to the NGC 6503 analysis. 

3.1 Kinematic Models 

The kinematic models implem ented in DiskFit are the same 
as i n velfit, as described by Spekkens fc Sellwoodl (|2007r i 
and ISellwood fc Sanchez! (|2QloF )7 with the additional capa- 
bility of fitting a symmetric outer disc warp. The current im- 
plementation of DiskFit can fit the following model includ- 
ing rotation and non-circu lar flows of fixed azimuthal p hase 
in a thin disc (c.f. eq. 5 in ISpekkens fc Sellwoodl ([20o3)): 



V m 



Vsys + sin i 



[A 



t cos f 



Vm,t cos(m9b) cos 9 — V m ,r sin(m#f,) sin^Jl) 



where Vt is the circular velocity, V m ,t and V m ,r are the tan- 
gential and radial components of non-circular flows with 
harmonic order m — 1 or m = 2 in the disc plane, 
9 and 9b are the azimuthal angles relative to the ma- 
jor axis and the non-circular flow axis, respectively, and 
i is the disc inclination. If m = 1, the model describes 
a lopsided flow; if m = 2 the model is bisymmetric, 
and describes a barred or elliptical flow. DiskFit there- 
fore fits disc-plane models to the data, taking care of 
the projection of model components of order m into sky- 
plane component s of order m' = m ± 1. Unlike tilted-ring 
techniques (e.g. ISchoenmakers. Franx. fc de Zeeuwl 1 19971 ; 



IWong. Blitz, fc Bosma 2004), the physical interpretation of 
the best-fitting DiskFit models is therefo re straightforward . 

For consistency with other work (e.g. lSimon et al]|2005l ; 
iTrachternach et al.ll2008h . the code can also fit radial flows, 
although such flows are physically poorly motivated. This 
model assumes m = distortions to the flow in t he disc 
plane, given by eq. 7 in ISpekkens fc Sellwoodl (^OOTf ): 



V sys + sini [ Vt cos 6 + V r sin#] , 



(->) 



where V r is the radial flow component. 

DiskFit requires a flat inner disc, but unlike its prede- 
cessor velfit, it allows for a symmetric warp in the outer 



disc. The disc is assumed flat out to warp radius r w , beyond 
which both the ellipticity and the position angle of the line 
of nodes vary in proportion to (r — r w ) 2 . DiskFit can find 
the best fitting inner warp radius r w and peak change in 
ellipticity and position angle of the line of nodes, or it can 
hold any combination of these parameters fixed. 

3.2 Photometric Models 

The photometric mo dels applied by Di skFit are identical 
to those described by I Reese et al.l |2007l '). DiskFit fits pho- 
tometric images with up to 3 components: a disc, bar and 
bulge. DiskFit assumes a flat, intrinsically round disc and 
a linear bar with a different apparent (projected) ellipticity 
and position angle in the plane of the disc. Both the disc 
and bar light profiles are non-parametric, and thus do not 
have a specific functional form. 

The bulge is parametrized by the Sersic function 



T(r) = To exp < -B n 



l/r, 



1 



(3) 



where r e is an effective radius, n in the Sersic index, and 
To is the intensity scale. The constant B n is a function 
of the Sersic index, and is defined by the implicit relation 
T(2n) = 2 7 (2n, B n ) (|Grahamll200lh . We allow the bulge to 
be spheroidal, with the disc plane being the plane of sym- 
metry. 



4 RESULTS OF KINEMATIC MODELING 

In this section, we present the best-fitting kinematic models 
to the Ha, CO and Hi velocity fields of NGC 6503. For each 
of the velocity fields in Fig. [T] we use DiskFit to fit for ro- 
tation only, as well as radial (m = 0), lopsided (m = 1) and 
bisymmetric (m = 2) non-circular motions. We also search 
for a symmetric warp signature in the Hi velocity field. We 
use DiskFit to determine the galaxy centre, inclination, sys- 
temic velocity, and position angle of the kinematic major 
axis. The minimum ring spacing is determined by the spa- 
tial resolution of the data and is increased to optimize the 
number of data points per ring, if necessary. 

We generate 1000 bootstrap realizations of each ve- 
locity field to determine uncertainties on the model pa- 
rameters: we assume correlated residuals over 4-5" to de- 
rive uncertainties on the irregularly sampled Ha data 
ISpekkens fc Sellwoodl 120071 . the only available approach 
in DiskFi t for this data type) , and the radial/re-scaling 
method of ISellwood fc Sanchezl (|2010l ) for the CO and Hi 
data. The parameter Aism is added in quadrature to the 
uncertainties in the emissio n line centroids during the fit 
l|Spekkens fc Sellwoodl 120071 ) . It is both an indicator of the 
turbulence in the disc and the relative size of the errors on 
the velocity field points. We set the value of Aism to produce 
reasonable x 2 statistics in the rotation-only model, and keep 
this parameter fixed for all other fits. The relative differences 
in x 2 f° r each dataset thus indicate the goodness of fit. We 
check that the model parameters and their uncertainties are 
unaffected by the particular value of Aism or the bootstrap 
technique that is used. 

We plot the derived rotation curves in Fig. [5] and the 
DiskFit rotation-only kinematic models and residuals in 



6 Kuzio de Naray et al. 



150 



100 



50 



• DensePak - Ha 
o SparsePak - Ha 

_J , , i , I , , , i L_ 

50 100 150 

r (arcsec) 




40 60 
r (arcsec) 



150 



100 



50 



200 



400 600 
r (arcsec) 



800 



Figure 5. Derived Ha (top), CO (middle), and Hi (bottom) ro- 
tation curves for NGC 6503. We find no evidence for coherent 
radial, lopsided, or bisymmetric non-circular motions in any of 
the datasets, nor a symmetric warp in the Hi data. The vertical 
dott ed line in the bottom panel indicates the edge of the optical 
disc dde Vaucouleurs et al.|[l99ll ). (A colour version of this figure 
is available in the online journal.) 



Fig. [6] We also discuss the derived constraints on non- 
circular motions in the disc (Fig. [7]). We compare the results 
for each dataset (Fig. [8} and to previous results in the lit- 
erature. Finally, we combine our multi-wavelength data and 
derive a master rotation curve for NGC 6503 (Fig. The 
best-fitting rotation-only model parameters for each dataset 
are listed in Table [1] 



4.1 Ha rotation curves 

We begin by fitting rotation-only models to both the 
DensePak and SparsePak velocity fields, and find that they 
provide a good representation of the data. The rotation ve- 
locities derived from the independent Ha datasets agree 
well with each other, particularly for 10" < r < 80" 
(see Fig. [5]l . The DensePak and SparsePak rotation curves 
both show a slow and smooth rise before flattening off at 
V ~ 120 kms -1 . The derived values of the position angle 
(DPak: -60.5 ± 0.6°; SPak: -59.6 ± 0.7°) and systemic ve- 
locity (DPak: 28.4 ± 1.0 kms -1 ; SPak: 28.4 ± 1.7 kms" 1 ) 



are identical to within their uncertainties. The best-fitting 
inclinations (DPak: 71.0±0.9°; SPak: 67.5±1.3°) however, 
are different at about the 2a level. Given that the radial 
extent of the SparsePak data is ~1.5 times smaller than the 
DensePak data, these differences are not unexpected. In ad- 
dition, there are approximately 3.5 times more data points 
in the DensePak velocity field than the SparsePak velocity 
field. That the adopted values of Aism for the two Ha ve- 
locity fields are different (6.5 kms -1 and 4.5 kms - for the 
DensePak and SparsePak data, respectively) is a reflection 
of the uncertainties on the individual fibre velocities being 
estimated in different ways for each dataset rather than an 
indication of a physical difference in the disc. 

Because there is such broad agreement between the two 
derived rotation curves and sets of model parameters, we 
combine the DensePak and SparsePak velocity field data. 
In Fig. [6] we show the combined Ha DiskFit rotation-only 
model and residuals. We also fit m = (0, 1, 2) models to the 
combined Ha velocity field; the tangential and radial com- 
ponents of these models are plotted in Fig. [7] Radial, lop- 
sided, and bisymmetric flows are consistent with kms -1 
throughout the optical disc, the change in \ 2 when these 
flows are included is not statistically significant, and the 
best-fitting disc parameters remain essentially unchanged 
from the rotation-only values in Table [T] This suggests that 
the best physical model for the Ha kinematics is one that 
con tains rotation only. 

lEpinat et all (|2008l ) observed NGC 6503 as part of the 
GHASP survey. Our final IFU Ha rotation curve agrees 
well with the Ha rotation curve they derived from their 
Fabry-Perot observations (see their figure E14). Both ro- 
tation curves extend to 150" and reach similar maximum 
velocities (IFU: 126±3kms -1 ; Fabry-Perot: 117±9kms -1 ). 
Their derived disc inclination (72 ±2°) and position angle 
(—61 ± 2°) are also both in very good agreement with the 
values determined from our IFU data (see Table [!}. 



4.2 CO rotation curve 

In our best-fitting rotation-only model of the CO velocity 
field, the CO rotation curve (see Fig. shows a smooth and 
steady rise in velocity (Vmax = 117 ± 8 kms -1 ) out to the 
last rotation curve point (Rmax = 80"). The best-fitting disc 
parameters for this model are given in Table [T] Fig. [7] shows 
the tangential and radial components of m = 0, 1, 2 models 
fit to the data. Similar to the Ha data, the CO kinematics 
lack signs of organized radial, lopsided, or bisymmetric flows, 
and we therefore find the rotation-only model to be the best 
des cription of the CO velocity fi e ld (se e Fig. [6} . 

iNishivama. Nakai. fc Kund (|200ll ) derive a CO rota- 
tion curve using 16" resolution observations from the 45-m 
Nobeyama Radio Observatory. Their data show rising ve- 
locities out to a radius of ~ 55", followed by a region of 
constant velocity of ~ 125 kms" 1 out to the end of the ro- 
tation curve just beyond 80". They determine the position 
angle of the major axis to be —57°, the inclination to be 
70°, and the maximum rotation velocity to be 127 kms -1 . 
Ou r CO results for NGC 6 503 are in good agreement with 
the INishivama et ail l|200ll ) data. 
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Figure 6. Rotation-only DiskFit models (top row) and residuals (bottom row) for the combined Ha (left column), CO (middle column), 
and Hi (right column) data. Note that the angular scale of each column is different: the horizontal bars in each panel of the top row 
are all 50" (1.2 kpc) long. See Fig. [T]for the explicit dimensions of each velocity field. (A colour version of this figure is available in the 
online journal.) 



4.3 Hi rotation curve 

The Hi rotation curve (see Fig. [5} rises to a maximum ve- 
locity of ~130 kms -1 around 100" and then slowly declines 
out to ~275" , just beyond the edge of the optical disc, be- 
fore flattening off with a rotation speed of ~115 kms - . 
Al though we rec o ver th e same rotation speed as reported 
by IGreisen et all (|2009l ). those authors do not find a grad- 
ual decline at 100" < r < 200": this difference stems from 
our adoption of an envelope tracing model to estimate the 
velocity field from the Hi data cube rather than the first mo- 
ment (Section 12. 3[1 . Beyond 600", the rotation curve shows 
hints of a very gradual decline , in rough agreement with the 
results of IGreisen et al.l (|2009h . 

We find radial, lopsided, and bisymmetric flows in the 
Hi data that are consistent with kms -1 (see Fig. [7}. The 
DiskFit velocity field model and residuals are shown in 
Fig. [6] like the optical and millimetre kinematic data, a 
rotation-only model provides the best description of the Hi 
velocity field. The sharp line along the minor axis in the 
residual map is an artefact introduced during the Modified 
Envelope Tracing (see Section I2.3[l . We also use DiskFit 
to search for a symmetric warp in the Hi disc, and find no 
evidence for one: specifically, we find a 3a upper limit on 
the maximum change in position angle at r > 30" to be 
wphim = 3°, while that on the disc ellipticity is welm =0.09. 
The gradual decline in the Hi rotation curve for r > 600" is 
therefore not caused by a warp. 

As in the models of IGreisen et all (|2009l ). our best- 
fitting rotation-only Hi model underestimates the projected 
velocity in the outermost regions in the south/southwest 
quadrant (bottom right) of the disc by 10 — 15 kms -1 (see 
Fig. |6]) . The position- velocity diagrams of IGreisen et al.l 



(|2009h suggest that this feature is a local warp in the disc. 
That DiskFit does not recover this feature as a warp is not 
surprising given its lack of symmetry and relatively small 
spatial extent. 



4.4 Comparison of velocity data 

The rotation curves for NGC 6503 derived using different 
tracers are all well-behaved and show no signs of radial (m = 
0), lopsided (m = 1) or bisymmetric (rn = 2) non-circular 
motions. In Figure [5] we directly compare the Ha, CO, and 
Hi rotation curves. The rising part of the rotation curve is 
well-mapped by the CO and Ha data. The overlap between 
the CO and Ha rotation velocities suggests that the optical 
data out to ~80" are n ot significantly affected by extinction 
fe.g. lBosma et al.lll992l ). Beyond ~45" , the Hi rotation curve 
agrees well with the Ha and CO data; the discrepant Hi 
rotation curve point at r ~ 30" is unreliable because it is 
likely biased by beam smearing. 

The position angle of the kinematic major axis is well- 
constrained, with all three datasets providing nearly in- 
distinguishable results (combined Ha: —60.3 ± 0.5°; CO: 
-60.5 ± 0.5°; Hi: -60.2 ± 0.3°). The disc inclinations de- 
rived from the CO (73.7 ± 1.1°) and Hi (73.5 ± 0.4°) are 
in good agreement, suggesting the inner and outer discs are 
coplanar. This is consistent with DiskFit not finding a warp 
in the Hi data (see Section 14. 3p . The Ha-derived inclination 
(70.0 ± 0.8°), however, is lower by about 2.5<j. There is a 
similar difference between the Hi (24.7 ± 0.7 kms -1 ) and 
Ha (28.5 ± 0.9 kms -1 ) systemic velocities. 

The agreement between the best-fitting model param- 
eters and similarities between the rotation curves for each 



8 Kuzio de Naray et al. 



Table 1. Best-fitting kinematic parameters for each dataset. Columns 2 and 3: position of centre relative to 17 h 49 m 26.4 s +70°08'40". 
Column 4: position angle of the kinematic major axis. Column 5: inclination assuming the thin disc approximation. The inclinations 
increase to a maximum of ~ 77° if an intrinsic axis ratio of 0.2 is assumed. Column 6: heliocentric systemic velocity. Column 7: ISM gas 
turbulence adopted in the fit. Column 8: Reduced chi-square of the best fitting model. The last row of the table gives the parameters 
adopted in deriving the multi- wavelength master rotation curve presented in Section 14.41 
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Figure 7. Top: Radial (m = 0) motions determined for the com- 
bined Ha (circles), CO (squares), and Hi (stars) data. Middle: 
Tangential (open points) and radial (filled points) components of 
the lopsided (m = 1) model. Bottom: Tangential (open points) 
and radial (filled points) components of the bisymmetric (m = 2) 
model. (A colour version of this figure is available in the online 
journal.) 
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Figure 8. Comparison of the derived rotation curves for 
NGC 6503 out to the edge of the optical disc. The vertical ar- 
rows at 3.5" and 37.5" denote the positions of the circumnuclear 
disc and star-forming ring, respectively, that are discussed in Sec- 
tion [5T2] (A colour version of this figure is available in the online 
journal.) 
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dataset suggest that all are tracing the same underlying disc 
kinematics. For these reasons, we determine a master rota- 
tion curve for NGC 6503 by combining the higher-resolution 
Ha data across the inner disc and the Hi data farther out; 
we use the Ha-derived galaxy centre and the Hi-derived in- 
clination, position angle, and systemic velocity. The master 
rotation curve is shown in Figure and the adopted disc 
geometry is given in Table [1] The rotation curve data are 
given in Table |A"T1 



Figure 9. Master rotation curve for NGC 6503. The inner rota- 
tion curve is determined from the Ha data (open circles) while 
the rotation at large radii is probed by Hi (stars). The rotation 
curve data are given in Table IATI The vertical dashed line shows 
the edge of the optical disc. 

5 RESULTS OF PHOTOMETRIC MODELING 

In this section, we present the best-fitting photometric mod- 
els to the Jfs-band image of NGC 6503. We use DiskFit 
to fit disc-only, disc+bulge, disc+bar, and disc+bar+bulge 
models. DiskFit determines the galaxy centre, disc incli- 
nation and position angle, the Sersic index, effective radius 
and ellipticity of the bulge, and the bar position angle and 




Figure 10. Masked, sky-subtracted _R"s-band image of NGC 6503 
(top), the best-fitting DiskFit disc+bar+bulge model (middle), 
and the residuals (bottom). Tickmarks in the top panel are spaced 
every 15" and the image has been rotated 30° to the west. The 
full radial extent of the CFHT image is ~ 370"; the panels are 
cropped to exclude most of the sky. The units on the colourbars 
are ADU. (A colour version of this figure is available in the online 
journal.) 

cllipticity. In addition, DiskFit calculates the percentage of 
light coming from each component (disc, bulge, bar) of the 
model. We generate 1000 bootstrap realizations of the image 
to determine uncertainties on the model parameters. 

We plot the best-fitting DiskFit models and residuals 
in Figures [TUll 111 and 1 121 the values of the model parameters 
are listed in Table [2] We discuss the photometric signatures 
of morphological features in the disc (Fig. I13[) . Finally, we 
present the derived Ks-band surface brightness profile for 
NGC 6503 in Figures [T4l and [TBI and compare it to previous 
results in the literature. 



5.1 DiskFit Photometric Models 

With different combinations of disc, bulge and bar com- 
ponents, we fit four DiskFit models to the CFHT data 
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Figure 11. Residuals for the disc-only (top), disc+bulgc (middle) 
and disc+bar (bottom) models. The size of each panel is the same 
as in Fig. 1101 the horizontal bar in the lower left corner of the top 
panel is 15" long and the models have been rotated 30° to the 
west. The units on the colourbar are ADU. (A colour version of 
this figure is available in the online journal.) 

(see Figures 1101 1111 and I12|l . The derived galaxy centre and 
the position angle of the disc are identical to within their 
uncertainties for all four models (see Table [2]). These disc 
position angles are also in very good agreement with the 
kinematically-derived position angle (see Table [JJ. With the 
exception of the disc-only model, the derived disc inclina- 
tions are, within their errors, consistent with each other and 
the kinematically-derived inclination. 

The disc-only model is a poor fit to the data. The de- 
rived disc inclination is comparatively low, and the residuals 
are large in the outer parts of the disc, as well as near the 
centre (see Fig. Illj l . The low inclination is an indication that 
there is a relatively bright circular component /feature at the 
galaxy centre that the model is trying to fit. In its innermost 
regions, NGC 6503 is more than a pure disc galaxy. 

The two-component (disc+bulge and disc+bar) models 
are an improvement over the disc-only model. In both of 
these cases, the disc contributes the majority of the light 
(~ 95%) and the models do a good, and similar, job of fit- 
ting the outer disc (see Fig. [TTJ . The smaller residuals in 
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Fig. [11] show this qualitatively, as do the smaller Xr values 
for the disc+bulge and disc+bar models in Table [2] the im- 
provement in the fit is highly statistically significant, even 
accounting for the extra parameters in the fit. 

With all of the bulge parameters (r e , n, e) in the 
disc+bulge model allowed to vary, the Sersic index of the 
bulge is poorly constrained. We therefore run two models 
and fix the bulge Sersic index to n = 1.1 an d n = 1.9. We 
choose these values based on the results of IPuglielli et al.l 
wh o model t he B- and i?-band surface brightness 
profiles of iBottemal (| 19891 ). They find the galaxy bulge to 
be nearly exponential (n = 1.1, r e = 6.4") or to be in the 
range of a pseudobulge [n = 1.9, r e — 14"), depending on 
whether they model the entire surface brightness profile, or 
exclude the Type II hump (see the discussion in Section |5, 31 
below). 

The results for the n = 1.1 and n = 1.9 DiskFit mod- 
els are both qualitatively and quantitatively very similar, 
and the DlSKFlT-derived effective radius of t he bulge when 
n = 1 .1 is the same as that determined by IPuglielli et all 
(|2010l ) to within the errors. We list the values of the n = 1.1 
model parameters in Tableland note here that, to within 
the errors, the disc parameters of the n = 1.9 model are 
the same as those derived in the n = 1.1 model and that 
the n = 1.9 bulge is slightly larger (r e = 8.9") and more 
round (e = 0.34) than the n = 1.1 bulge. For concreteness, 
we adopt the n = 1.1 bulge model parameters in what fol- 
lows, although the analysis remains unchanged if the n = 1.9 
results are used. 

In Fig. \12k we show the n = 1.1 DiskFit bulge. In the 
disc+bulge model, the bulge is not spherical, but rather has 
a projected ellipticity along the disc major axis of e — 0.4. It 
contains 5.4% of the total light and has an effective radius 
of r e = 6.7". 

The bar determined by DiskFit in the disc+bar model 
is relatively round (e = 0.39) and oriented ~ 7° away from 
the disc major axis in projection. It contains 4.5% of the 
total light. The bar is shown in Fig. 112b . 

The bar model does a better job of fitting the data 
than the bulge model {xl,bar = 2 - 72 vs - xl,bui ae = 3.05). 
This can also be seen by looking closely at the residuals in 
Fig. I12H and Fig. I12fe . In particular, in the bulge model, 
the large ±10 ADU residuals at the centre, the red/pink 
regions of 5-8 ADU residuals in the upper right quadrant 
and the extended red region of 5 ADU residuals at the left 
side of the panel of Fig. I12H are greatly improved in the bar 
model (Fig. 112b ). The trade-off in the bar model are larger 
residuals 90° away from the bar in the upper left and lower 
right quadrants. 

The similarity between the morphology of the bulge in 
the disc+bulge model and that of the bar in the disc+bar 
model (Fig. 112b and Fig. 112b ) illustrates that the bar com- 
ponent is trying to encompass the bulge of the galaxy, while 
simultaneously trying to fit an asymmetric component that 
is not aligned with the major axis of the disc. The non-disc 
components of the disc+bulge model and disc+bar model 
are therefore both attempting to fit the same feature in the 
light distribution (Fig. I12[) . The primary difference between 
the models is whether or not the asymmetry is allowed to 
move away from the disc major axis. 

We are able to improve the fit in the entire central re- 
gion of the galaxy by allowing the DiskFit model to con- 
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Figure 12. (a): bulge component of the disc+bulge model, (b): 
bar component of the disc+bar model, (c): bar and bulge com- 
ponents of the disc+bar+bulgc model. The bar is the structure 
stretching from the lower left corner to the upper right, (d): 
residuals of the disc+bulge model, (e): residuals of the disc+bar 
model, (f): residuals of the disc+bar+bulge model. Each panel is 
25" x25" in size and the models have been rotated 30° to the west 
so that the disc major axis is horizontal. The tickmarks in panel a 
are spaced every 5". The units on the upper and lower colourbars 
are ADU. In the two-component models (panels a and b), the 
non-disc structures are very similar except for their position an- 
gles relative to the disc major axis. In the three-component model 
(panel c), the bar is more elliptical and angled farther away from 
the disc major axis than the bar in the disc+bar model (panel b). 



tain both a bar and a bulge (see Fig. [I0]and Fig. [12}. In this 
three-component model, the disc contains 94.7% of the total 
light and its geometry agrees well with those derived from 
the kinematic data. The bar contains 0.7% of the light, has 
an ellipticity of e = 0.53, and projects near the minor axis of 
the disc (PAf, ar ~ —17°). As the bulge parameters are un- 
constrained if they are all allowed to vary, we fix the values 
to those found in the n = 1.1 disc+bulge model. Changing 
the bulge properties does not make a difference to the Disk- 
Fit model: it is the addition of an extra physical component 
that matters. The bulge contains 4.6% of the total light. 

The modest improvement in Xr f° r the three-component 
disc+bar+bulge model (see Table [2} is statistically signifi- 
cant at the 4ct level due to the large number (> 37500) 
of degrees of freedom. Moreover, we consider the three- 
component model to be the best physical description of 
NGC 6503 because the residuals in Fig. [^highlight the need 
for both a component in the disc and an extended asymmet- 
ric structure that is not aligned with the disc major axis. 

5.2 Photometric Signatures of Morphological 
Features 

Using a multi- wavelength dataset, iFlOl determine that the 
disc of NGC 6503 contains several structures. They find that 
there is a 3.5" circumnuclear disc that is closely aligned to 
the major axis of the galaxy disc and that a nuclear spiral 
exists within the circumnuclear disc. They also find a 7" 
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Figure 13. Observed if s-band image of NGC 6503 (top) and the 
individual components of the DiskFit disc+bar+bulge model. 
The large and small dashed ellipses indicate the positions of the 
star-forming ring and circumnuclear disc, respectively, that are 
discussed bv lFlOl , To guide the eye, the dashed lines are parallel 
to the spiral arms coming off the ends of the end-on bar; for 
comparison to Fig. 1101 these lines are 15" (0.35 kpc) in length. 
The image has been rotated 30° to the west. (A colour version of 
this figure is available in the online journal.) 



end-on bar; the galaxy's spiral arms extend from the bar 
ends. Also present is a 37.5" star-forming ring. 

In Fig. 1131 we indicate the positions of the features 
highlighted by Flfj in relation to the disc, bar and bulge 
components determined by DiskFit in the three-component 
model. The DiskFit bulge is clearly fitting the circumnu- 
clear disc. DiskFit finds a relatively bright bar-like compo- 
nent within the circumnuclear disc, although the negative 
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Figure 14. Surface brightness profile for NGC 6503 derived from 
the DiskFit disc+bar+bulge model. The solid black line is the 
total light. The blue open circles are the contribution from the 
disc. The red open triangles are the contribution from the bar. 
The green stars are the contribution from the bulge. The arrow 
at 37.5" indicates the position of the star-forming ring. (A colour 
version of this figure is available in the online journal.) 



residuals in this region (see Fig. I12F ) suggest this to be an 
artefact of the fit. On the other hand, the faint bar-like com- 
ponent beyond the circumnuclear disc encloses an elliptical 
area that roughly connects the galaxy's spiral arms: this 
recovered structure is likely real. We note that Fig. [13] has 
been rotated 30° to the west making the close-to-end-on bar 
(PAbar ~ —17°) appear closer to ~ —45° in that figure. The 
star-forming ring is not prominent in the DiskFit model 
nor the Tvs-band image, which is not surprising given that 
the Tvs-band probes the old stellar population. 

The ivs-band photometric modeling with DiskFit 
shows NGC 6503 to be a disc-dominated system that has 
small, but measurable, non-disc structure. To summarize, 
the disc+bar+bulge model described in Section 15.11 is not 
only the best statistical description of the data, but also 
th e bes t physical model when compared to structures found 
bv lFldl . To accurately model the inner region of the galaxy 
both on and off the disc major axis, two separate non-disc 
components are required. In addition, the bar that D iskFi t 
finds is positioned very similarly to the bar found bv lFlOl . 

5.3 Surface Brightness Profile 

In Fig. 1141 we show the i\s-band surface brightness pro- 
files derived from this disc+bar+bulge model. Note that 
DiskFit automatically produces separate profiles for each 
of the model components. The high central point of the bar 
profile is unphysical and is likely being confused with the 
bulge component. This means that the fraction of the to- 
tal light that is assigned to the bar (0.7%) is actually an 
over-estimate. 

The radial extent of the surface brightness profile in 
Fi g. 1141 is three times farther than previous data compiled 
bv iFlCl. The pr ofile displays a Freeman Type II behaviour 
|Freemanlll970l ). being exponential at large radii, having a 
plateau between 10" and 40", and rising steeply inside of 
10". To highlight the plateau, in Fig. [15] we fit an expo- 
nential to the profile at radii between ~ 40" and 120" and 
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Table 2. Best-fitting photometric parameters. Rows 1 and 2: position of centre relative to 17 h 49 m 26.4 s +70°08'40". Row 3: position 
angle of the disc major axis. Row 4: inclination of the disc assuming the thin disc approximation. The inclinations increase to a maximum 
of ~ 77° if an intrinsic axis ratio of 0.2 is assumed. Row 5: effective radius of the bulge. Row 6: Adopted Sersic index of the bulge. Row 
7: ellipticity of the bulge. Row 8: position angle of the bar in the sky plane. Row 9: ellipticity, defined as l-(b/a), of the bar. Row 10: 
reduced chi-square of the model. Row 11: fraction of light attributed to the disc. Row 12: fraction of light attributed to the bulge. Row 
13: fraction of light attributed to the bar. 
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Figure 15. Top panel: total inner DiskFit surface brightness 
profile for NGC 6503. The red line is an exponential fit to the 
data between 40" and 120". Bottom panel: difference between 
the exponential fit and the observed surface brightness profile. 
The vertic al arrows indicate the position of the star-forming ring 
from iFld that coincides with the Type II plateau in the surface 
brightness profile. (A colour version of this figure is available in 
the online journal.) 



plot the difference between the exponential fit and the ob- 
served profile. There is a maximum difference between the 
two of about 0.8 mag arcsec -2 at a radius of ~ 10". The 
position of the star- forming ring discussed bv lFld coincides 
with the outer radius of the plateau in the surface bright- 
ness profile. The outer profile beyond 120" (Fig. 1 14 p is not 
exponential. We have verified that this is a real feature of 



the galaxy and not an artefact related to sky subtraction 
(e.g. lMcDonald. Courteau. fc Tullvll2009t ) . 

present an iiT-band surface brightness profile of 
NGC 6503, along with B- and i?-band profiles from lBottemal 
l|l989f ). and 3.6/xm and 4.5/xm profiles from the Spitzer 
archive. These profiles show the same Type II behaviour 
at the same radii as our Tfs-band profile. In addition, the 
magnitude of the average deviation of these data from an 
exponential profile is consistent with what we find for the 
Jfs-band data. 



DOES NGC 6503 CONTAIN A STRONG 
BAR? 



iBottemal (|l989T l argues that dust extinction is the cause of 
the plateau observed in the B- and J?-band surface bright- 
ness profiles of NGC 6503. We can rule out this explana- 
tion as not only is the Type II profile observed in our in- 
frared 7\s-band data, but the agreement between our de- 
rived Ha and CO rotation curves a lso suggest s minimal ex- 
tinction. Based on simulations by lB05l . IFld argue that a 
better explanation for the Type II sur face brightness profile 
in NGC 6503 is the presence of a bar. IPuglielli et all (|2010h 
similarly prefer th e bar explanation based on their m odel- 
ing and results of Fovle. Courteau. fc Thackerl (|200ct ) and 
iBournaud. Combes, fc Semelinl ( 20051 ). 

IFld provide a number of lines of evidence to support 
the idea of a bar. Nuclear spirals and star-forming rings 
are r esonance features of bars (e.g. lAthanassoula fc Bureau! 
Il999h . and both are observed in NGC 6503. There is 
a significant central drop (cr-dro p) in the stel lar veloc- 
ity dispersion of the galaxy (e.g. IBottemal 1 19891 ). cr-drops 
can be caused by a number of mechanisms including 
bars, dynamically cold gas forming stars in a circu mnu- 
clear disc, or even counter-rotating s tellar di scs (e .g. lB05l : 
lEmsellem et al. 200 ll : iRix et al. 1991 ). While [Fig find ev- 
idence for a circumnuclear disc in their WHIRC ff-band 
imaging data, they argue that the cr-drop is more likely 
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caused by a bar, given the presence of the nuclear spiral 
and star- forming ring. 

iFlOl specifically indicate that it is a strong end-on 
bar. |B05| run TV-body simulations to highlight how the bar 
strength and orientation (with respect to the observer's line- 
of-sight) and galaxy inclination affect, among other things, 
the shape of the a profile. In these simulation s, onl y a strong 
end-on bar can cause a central cr-drop. As in lFlOl . we adopt 
this qualitative distinction between a strong bar and a weak 
one for the following discussion. 

Do we find any direct evidence for a strong bar in our 
dataset and DiskFit models? Strong bars should pr oduce 
disturbed kinematics similar to those seen in figure 1 of lB05l . 
Our Ha, CO, and Hi velocity fields and derived rotation 
curves are well-described by rotation-only DiskFit models 
and agree remarkably well with each other. With the ex- 
ception of a possible hint of the broadening of the velocity 
profile in the central pixel of the CO data, we do not detect 
the bar in the rotation curves or velocity fields. 

Though we can see hints of the star-forming ring in the 
Ha and CO intensity maps (see Figures [5] and [3J , we do 
not detect it kinematically (see Figures [6] and suggesting 
that the bar is not currently strong enough to be funneling 
gas. Additionally, there is no build-up of CO along the ring 
near the ends of the bar (see Fi g. EH), which one would ex - 
pect if the bar were strong (e.g. Lin. Yuan, fc Butall2008l ). 
We estimate the molecular gas mass at 30" < r < 40" 
— the location of the star-forming ring — to be 
4.2 x 10 6 M Q (using X C o = 2 x 10 20 cm~ 2 /K kms" 1 an d 
a C o = 4.3 Mq/K kms" 1 ; ISoIomon fc Vanden Boutl [20051 ). 
much less than than the enclosed mass of the galaxy at 
r~40", ~ 2.2 x I0 9 M Q . That we see a star-forming ring 
is not conclusive evidence for a bar; rings can be formed in 
a variety of ways including during accreti on events or col- 
lisions/interactions between galaxies (e.g., iBerentzen et al.l 
l2003l ; lMapelli fc Mave"3l2012h . By themselves, our kinematic 
data do not suggest that a bar is currently present in 
NGC 6503. 

While our data show no evidence for a bar in the kine- 
matics (beyond 5"), a faint, but statistically significant bar- 
like feature is detectable in the photometry. To fully model 
the light distribution in NGC 6503, both a bar and a bulge 
are required. Can this photometric feature consistently ex- 
plain the a-drop, Type II profile, and gas morphology given 
the regular kinematics? 

Given the faintness of the photometric bar and the very 
well-behaved kinematic data, we find it is most likely that 
this bar is weak and it is the combination of the bar and 
the circumnuclear disc in NGC 6503 that work together to 
produce the cr-drop. Without more stellar kinematic data, 
however, we cannot conclusively rule out the possibility that 
other mechanisms, as described above, are not responsible 
for creating the cr-drop. 

On the contrary, it is plausible that the observed Type 
II surface brightness profile could be caused by th e photo- 
metric bar. Consistent with the simulations of |B05| . the bar 
we detect could be strong enough to rearrange the mate- 
rial in the disc, producing t he ch aracteristic Type II profile. 
The i = 75° simulations of |B05| (see their figure 4) and an 
intrinsic axial ratio q ~ 0.2 (Tables Q] and [S]), appropriate for 
the disc geometry of NGC 6503, suggest that relatively weak 
bars may produce this feature because the low-density minor 



axis of the bar is exposed. The bar also has to be very close 
to end-on (as the DiskFit model suggests it is), otherwise 
the characteristic signatures in the surface brightness profile 
would be less pronounced and the galaxy isophotes would 
have more of a "peanut" -shape rather than a "bulge+ring"- 
shape (see the right-most column of figure 4 in lB05l ). 

It may be possible to explain the properties of 
NGC 6503 by invoking a strong bar which dissolved some- 
time in the past. For this scenario to work in NGC 6503, the 
bar must have once been much stronger, so as to be able to 
rearrange the disc material, before significantly weakening, 
or nearly dissolving, so that we now see very little trac e of it 
in either the kinematics or photometry. According to |B05| . 
who also study the time evolution of the surface brightness 
and stellar velocity dispersion profiles, it is only late in the 
lifetime of the bar, well after it has formed and buckled, 
that the Type II plateau and central cr-drop become signif- 
icant and well-established. This would suggest that the bar 
in NGC 6503 is "old" . It also suggests that the star-forming 
ring is not caused by the bar (see above for alternative ring 
formation mechanisms), as it would be puzzling for such a 
bar to now be re spon sible for the formation of a young (less 
than ~ 0.5 Gyr; IFlOl ) ring. Estimating the decay timescale 
for a bar is a challenging numerical problem and different 
techniques and studies yield conflicting results, though the 
general co nclusion is that ba r s can weaken and/or disappear 
(see, e.g., IShen fc Sellwoodl |2004| ; iKormendv fc Kennicuttl 
|2004 for a discussion). 

In the absence of a full suite of models, we cannot con- 
clusively determine the evolutionary history of NGC 6503 
and unequivocally establish that a (strong) bar is respon- 
sible for the observed cr-drop, Type II surface brightness 
profile and star-forming ring. We can say, though, that a 
strong bar is not currently present in the galaxy. The bar 
that exists today in NGC 6503 is weak enough to avoid be- 
ing detected in the gas kinematics, but is just strong enough 
to leave a measurable photometric signature. 



7 SUMMARY 

We have obtained new, multi- wavelength, high-quality spec- 
troscopic and photometric data for the nearby spiral galaxy 
NGC 6503. Specifically, we have presented two Ha veloc- 
ity fields obtained using the DensePak and SparsePak IFUs, 
a CO(J=I-0) velocity field obtained using CARMA, and a 
re-derived VLA Hi velocity field. We have also presented a 
CFHT WIRCAM A's-band image of the galaxy. 

We have introduced the first public release of Disk- 
Fit, a code that fits non-parametric models to either veloc- 
ity fields or images. The kinematic models implemented in 
DiskFit fit for rotation, as well as m — 0, 1, 2 non-circular 
flows. DiskFit is also able to fit for a symmetric outer disc 
warp. In the photometric models, DiskFit can fit up to 3 
components: a non-parametric disc and bar, as well as a 
Sersic bulge. Relative to other algorithms, DiskFit excels 
at finding relatively weak, coherent asymmetries and pro- 
viding a robust estimate of the underlying disc properties in 
their presence. 

We have used DiskFit to model our velocity field data 
and search for non-circular motions. We find the velocity 
fields and derived rotation curves to all be well-described by 
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rotation-only models. We find no evidence for coherent ra- 
dial, lopsided, or bisymmetric non-circular motions in any of 
the datasets, nor a symmetric warp in the Hi data. We com- 
pare the best-fitting model parameters and rotation curves 
for each dataset and find them to be in good agreement 
with each other, suggesting that they are all tracing the 
same underlying disc kinematics. We combine the Ha and 
Hi rotation curves to determine a master rotation curve for 
NGC 6503. 

We have also used DiskFit to analyze the Tfs-band 
image of NGC 6503 and produce models of the disc, bar, 
and bulge. We find the galaxy to be an axisymmetric, disc- 
dominated system that contains a bulge and weak end-on 
bar. Together, the bulge and bar contribute only 5.3% of 
the total light. We find the galaxy to have a Type II surface 
brightness profile that could have been caused by the bar. 

NGC 6503 is a disc galaxy with regular, well-behaved 
kinematics. Despite indirect photometric indications of the 
presence of a bar (e.g., star-forming ring, Type II sur- 
face brightness profile) and the faint (0.7% of the total 
light) detection in our ifs-band photometry, we do not find 
kinematic counterparts to these features in the DiskFit 
models of the velocity field data; at least at present, the 
bar is weak. The large radial extent and gently declining 
outer rotation curve of the galaxy can place interesting con- 
straints on the concentration of the dark matter halo (e.g. 
ICasertano fe van Gorkom|[l99ll ); NGC 6503 thus remains 
one of the best targets for detailed studies of the dark and 
baryonic structure of late-type disc galaxies. However, par- 
ticularly in the event that the bar was once stronger and 
caused the Type II profile, to constrain the inner dark mat- 
ter halo, full modeling of a barred galaxy that follows the 
angular momentum exchange between the disc and halo is 
required. 
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APPENDIX A: MASTER ROTATION CURVE 

In Table I All we present the data for the master rotation 
curve for NGC 6503 discussed in Section l4~4l and plotted in 
Figure [5] 



Table Al. Master Rotation Curve for NGC 6503. 
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